Convective Conditions on MHD Casson Fluid in a Suspension of Dusty Particles with Darcy’s Porous

 

S. U. Mamatha1, C. S. K. Raju2*, D Madhavi Giridhar 3

1Department of Mathematics, Garden City University, Bengaluru, Karnataka 560049

2Department of Mathematics, GITAM University, Bangalore, Karnataka-562163

3Department of Mathematics, Garden City University, Bengaluru, Karnataka 577002

*Corresponding Author E-mail:mamathasupadhya@gmail.com1, sivaphd90@gmail.com2*, madhavi.g@gardencitycollege.edu

 

ABSTRACT:

The present study is to analyze Darcy-Forchheimer steady incompressible (MHD) non-Newtonian Casson fluid with suspension of dust particles over a stretching surface with exponentially decaying heat source. Convective boundary conditions are also considered. Similarity transformations were used to convert partial differential equations (PDEs) to a system of nonlinear ordinary differential equations (NODEs) which are solved numerically by employing Runge-Kutta with shooting technique. The effect of pertinent parameters on velocity and temperature profiles of both fluid and dust phase within the boundary layer has been studied by considering various values of controlling parameters. Additionally, skin friction coefficient and Nusselt number have been examined for various values of the governing parameters. It is found that non-uniform inertia coefficient of porous medium, porosity parameter and Casson fluid parameter decreases friction factor and Nusselt number.

 

KEYWORDS: Exponentially decaying heat source, Convective conditions, Dusty fluid, Casson fluid, Darcy-Forchheimer.

 


INTRODUCTION:

Over recent years, researchers are allured by the characteristics and performance of non-Newtonian fluids. Such motivation is due to complex physical behavior of non-Newtonian fluids. A non-Newtonian fluid finds numerous applications in pharmaceuticals, food products, crystal growth, fiber technology, chemical and material engineering. Hence several constitutive models for non-Newtonian fluids are available in the existing literature [1-4]. Non- Newtonian Casson fluid is a shear thinning fluid which exhibit yield shear stress. If suppose the yield stress is lesser than the applied shear stress than the liquid would start moving. For example tomato sauce, honey, human blood and fruit juices.

 

Casson fluid finds several applications in cancer homeo-therapy, fibrinogen, protein cells and red blood cells.  Because of these applications researchers Raju et al.[5] studied heat and mass transfer of MHD Casson fluid on an exponential permeable stretching surface. Further  Raju et al. [6] took up a comparative study of non-Newtonian Casson and Carreau fluid over a stretching surface and found that heat and mass transfer in Casson fluid is considerable high compared with Carreau fluid. 

 

Phenomenon of suspended particles finds significant role in advanced processes of engineering and scientific field concerned with polymer technology, metallurgy, powder technology, power generation, atmospheric fallout, nuclear reactor cooling, combustion, waste water treatment, electrostatic precipitation of dust ,fluidization etc. Fundamental studies of suspended particles in the flow were initiated by Saffman [7]. Recently Krupa Lakshmi et al.[8] took up heat and mass transfer of a dusty liquid past a stretching sheet  by considering thermal radiation in the energy equation. Gireesha et al. [8] studied influence of viscous dissipation on flow of dusty fluid over stretching sheet. Krishnamurthy et al. [9] took up theoretical study of suspended particle on slip flow and nanofluid in a porous Medium with the presence of nonlinear thermal radiation.

 

Extensive research has been undertaken by the researchers to study flow induced by a stretching surface in view of its large range of applications in chemical processing equipment, polymer extrusion, wire and fiber coating, glass fiber production, water pipes, irrigation channels etc. Due to its vast applications, initially, Sakiadis [11] presented the concept of two-dimensional boundary layer flows past a continuously stretching sheet with constant speed. Later on, several authors [12-22] extended the research over a stretching sheet for various fluid flow models over several geometries such as cone, Cylinder, Plate etc.

 

Heat and mass transportation from varied geometries embedded with porous media has several engineering and geophysical applications in geothermal reservoirs, thermal insulation, hydrology, agriculture, mineral processing, etc.  The classical Darcy’s theory which has been used in most of the flow problems is suitable for lower velocity and smaller porosity. When inertial and boundary effects take place with higher flow rate Darcy’s theory is inadequate. Keeping this in view researchers [23-25] studied heat and mass transfer of fluid by considering Darcy- Forchheimer porous medium.

 

By studying the above literature, in this paper we proposed a mathematical model of boundary layer flow of Casson fluid with suspended uniform sized dust particles towards a stretching sheet in the presence of magnetic field, Darcy-Forchheimer porous medium, exponentially decaying space dependent internal heat source and convective boundary condition. The governing nonlinear partial differential equations are reduced into a set of non-linear ordinary differential equations by using local similarity transformation. These nonlinear ordinary differential equations then solved numerically by Runge-Kutta based shooting process for different values of parameters of interest. Influence of governing parameters on velocity, temperature, Skin friction coefficient and local Nusselt number are analyzed and discussed in detail with help of plotted graphs and tables.


 

 

Fig.1 Schematic representation of the flow diagram

2. MODELING:

A steady two-dimensional, incompressible, magnetohydrodynamic boundary layer flow of Casson fluid with suspended uniform size dust particles over a stretching sheet characterizing Darcy-Forchheimerporous medium is developed. The stretching sheet is aligned with the   axis and the flow is confined to. The surface is stretched with the velocity  in direction with rate.Magnetic field of strength is imposed normal to the fluid flow as shown in Figure 1 By assuming electrical conductivity of the fluid to be small we have neglected the induced magnetic field. Convective surface temperature is characterized by and heat transfer coefficient by .The temperature  at the surface of the sheet is considered to be more than the ambient fluid temperature. The associated governing boundary layer equations are:

                                                                                                                                                         (1)

                                                  (2)

                                                                                                                                                      (3)

                                                                                                                              (4)

                 (5)

                                                                                                                 (6)

With boundary conditions 

                                                                                (7)

Here - velocity components along  direction of the fluid and dust particle phase, - velocity components along direction of the fluid and dust particle phase,are density of the fluid, density of the dust particle phase, mass of the dust particles, number density of the dust particle, kinematic viscosity, thermal conductivity, electric conductivity, uniform magnetic field, specific heat of the fluid and dust phase.-permeability of porous medium, A is power law index, denotes stokes drag coefficient. -radius of the dust particle. is the non-uniform inertia coefficient of porous medium, where is the drag coefficient [24] .represents temperature of the fluid, temperature of the dust particle phase, ambient fluid temperature, convective fluid temperature, and convective heat transfer coefficient. andrepresents relaxation time of dust particles and thermal equilibrium time.Coefficient of the dimensionless space-dependent internal heat generation [26]. 

Following similarity transformation were used to convert PDEs into set of ODEs.

                                                                              (8)

Continuity equations (1) and (3) are identically satisfied. Equations (2), (4)-(7) are transformed as follows.

   (9)

                                                                   (10)

      (11)

                                                                                                     (12)

The boundary conditions defined in equation (7) will be transformed to:

                                                                          (13)

 in the above expressions, represents magnetic parameter, mass concentration of dust particles, fluid-particle interaction parameter,Casson fluid parameter, porosity parameter, inertia coefficient, Prandtl number, heat source parameter, Eckert number, specific heat ratio, fluid particle interaction parameters,biot number.

                                                 (14)

Distinctive measures of practical interest are skin friction coefficient   and local Nusselt number  Which are defined as ;

Where  is the local Reynold’snumber.

 

 

 

 

 

Fig.2 The effect of on velocity field                                                    Fig.3 The effect of on temperature field

 

 

Fig.4 The effect oon temperature field                                                     Fig.5 The effect oon velocity field

 

Fig.6 The effect oon velocity field        Fig.7 The effect oon velocity field

 

 

Fig.8 The effect oon temperature field       Fig.9 The effect of  on temperature field


 

 

Fig.10 The effect o on temperature field

 

Table 1 Comparison of the results for local Nusselt number with.

Pr

Ishak

[17 ]

El-Aziz

[8 ]

K.L.Krupa Lakshmi

et.al[2]

Present Study

0.72

0.8086

0.80873

0.808630

0.8075

1.0

1.0000

1.0000

1.0000

1.0000

3.0

1.9237

1.92368

1.92367

1.9146

10.0

3.7207

3.7207

3.72067

3.6516

100

12.2941

12.2941

12.294087

12.2936

 

 

 

 

 

 

 


 

Table 2 Numerical values of Skin friction and Nusselt number for different values of .

0.5

 

 

 

 

 

-0.442845

-0.344656

1

 

 

 

 

 

-0.521758

-0.350769

1.5

 

 

 

 

 

-0.591849

-0.356023

 

0.1

 

 

 

 

-0.442845

-0.209401

 

0.3

 

 

 

 

-0.442845

-0.440393

 

0.5

 

 

 

 

-0.442848

-0.566926

 

 

0.1

 

 

 

-0.430626

-0.343955

 

 

0.3

 

 

 

-0.454752

-0.345334

 

 

0.5

 

 

 

-0.477702

-0.346621

 

 

 

0.1

 

 

-0.425760

-0.343328

 

 

 

0.3

 

 

-0.459462

-0.345950

 

 

 

0.5

 

 

-0.491397

-0.348431

 

 

 

 

0.1

 

-0.331109

-0.338235

 

 

 

 

0.3

 

-0.517415

-0.348701

 

 

 

 

0.5

 

-0.616068

-0.353715

 

 

 

 

 

0.1

-0.442845

-0.355359

 

 

 

 

 

0.3

-0.442845

-0.335014

 

 

 

 

 

0.5

-0.442845

-0.318266

 


RESULTS AND DISCUSSION:

For analyzing the approximate solutions of velocity as well as temperature fields, the non-linear ordinary differential Eqs. (9)-(12) with reference to the boundary conditions (13) are solved numerically with the assistance of Runge-Kutta and Newton’s methods. For numerical solutions we have considered the values of non- dimensional parameters asthese values are taken constant in this study besides the varied parameters as mentioned in the respective figure. In Figure 2and 3 the influence of magnetic field parameter on velocity and temperature profiles of fluid and dust phases can be observed. As accelerates, The Lorentz force which opposes the flow regime also improves hence there is  decrement in the and Lorentz force generates friction on the flow, which results more heat energy causing increase in the temperature distribution of the flow.

 

Figure 4 demonstrates the influence of convective heating parameter on the temperature distribution of fluid and dust phase. It is noticed that increasing values of improves thermal boundary layer. Physically convective heat exchange at the surface of the sheet helps to improve the thermal boundary layer thickness. Figure 5 presents larger values of non –uniform inertia coefficient  declines the velocity profiles of both fluid and dust phase.  Figure 6 and 7 depicts the variation in velocity and temperature profiles for larger values of porosity parameter. It is observed that increment in  reduces velocity field and improves temperature field for both fluid and dust phase. In Figure 8 and 9 it is noticed that increment in Casson fluid parameter decreases velocity boundary layer and improves thermal boundary layer for both fluid and dust phase. Figure 10 demonstrates the influence of exponentially decaying internal heat generation parameter  over temperature profiles of fluid and dust phase. It is observed that increase in  enhances heat generation which causes the temperature profiles of fluid and dust phase to enhance. 

 

Table 1 Shows the comparative study of thepresent result with the published results. Present results agree with the published results.

Table 2 display the effect of non-dimensional parameters on skin friction coefficient and local Nusselt number. It is noticed that rise in magnetic field parameter minimize the friction coefficient along with Nusselt number. Increment in Biot number decrease Nusselt number. But no variation is observed in friction factor. Rising values of inertia coefficient, porosity parameter and Casson fluid parameter decreases skin friction and Nusselt number. Increment in exponentially decaying internal heat generation parameter enhances heat transfer rate and friction factor remains same.

 

4. CONCLUSION:

Darcy –Forchheimer two- dimensional flow of non-Newtonian Casson fluid with suspended uniform size dust particles over a stretching sheet with exponentially decaying heat generation and convective boundary condition has been discussed. The major findings of present study are as follows.

·        The magnetic interaction parameter motivates friction factor coefficients and heat transfer rate.

·        Increment in Biot number improves the thermal boundary layer of both phases. But retards heat transfer rate.

·        Heat source parameter declines thermal boundary layer but improves heat transfer rate.

·        Non-uniform inertia coefficient, Porosity parameter, Casson fluid parameter reduces friction factor and heat transfer rate but improves thermal boundary layer in both the fluid and dust phase

 

5. REFERENCES:

1.     Raju CSK, Sandeep N, Jayachandra Babu M, Sugunamma V. Dual solutions for three-dimensional MHD flow of a nanofluid over a nonlinearly permeable stretching sheet. Alexandria Eng. J. 2016; 55:151–162. doi:10.1016/j.aej.2015.12.017.

2.     Mamatha SU, Mahesha, Raju CSK.Multiple Slips on drodynamic Carreau Dustynano Fluid Over a Stretched Surface with Cattaneo-Christov Heat Flux. J. Nanofluids. 2017;6: 1074–1081

3.     Mamatha S U, Raju CSK, Madhavi G , Mahesha. Unsteady 3D  MHD Carreau and Casson Fluids over a Stretching Sheet with  Non-Uniform Heat Source/Sink.IISTE. 2017;52: 10-23

4.     Raju CSK, Sandeep N, Sulochana N, Babu M.J. Dual Solutions of MHD Boundary Layer Flow past an Exponentially Stretching Sheet with Non-Uniform Heat Source / Sink.J. Appl. Fluid Mech. 2016; 9(2). 555-563 .doi:10.18869/acadpub.jafm.68.225.24784.

5.     Raju CSK, Sandeep N, Sugunamma V, Babu M.J, Reddy JVR,  Heat and mass transfer in magnetohydrodynamic Casson fluid over an exponentially permeable stretching surface. Eng. Sci. Technol. an Int. J. 2016 ;19 :45–52. doi:10.1016/j.jestch.2015.05.010.

6.     Raju CSK, Sandeep N, Unsteady three-dimensional flow of Casson –Carreau fluids past a stretching surface.AlexandriaEngineeringJournal.2016;55(2):1115-1126 http://dx.doi.org/10.1016/j.aej.2016.03.023.

7.     Saffman PG. On the stability of a laminar flow of a dusty gas. J. Fluid Mechanics.1962;131:20–128.

8.     Krupa Lakshmi KL, Gireesha BJ, Rama S R Gorla, Mahantesh B.Two-phase Boundary Layer Flow , Heat and Mass Transfer of a Dusty Liquid past a Stretching Sheet with Thermal Radiation. Int J. Industrial Mathematics .2016;8(3): 14 pages.

9.     Gireesha BJ, Roopa PG,.Bagewadi CS.Effect of  viscous dissipation and heat source on flow and heat transfer of dusty fluid over unsteady stretching sheet. Appl Math Mech. 2012; 33: 1001-1014.

10.   Krishnamurthy MR., Gireesha BJ, Gorla RSR, Prasannakumara B.C. Suspended Particle Effect on Slip Flow and Melting Heat Transfer of Nanofluid Over a Stretching Sheet Embedded in a Porous Medium in the Presence of Nonlinear Thermal Radiation. J. Nanofluids. 2016; 5:502–510.

11.   Sakiadis BC, Boundary Layer Behavior on Continuous Solid Surfaces:II Boundary Layyer on a Continuous Flat Surface. AIChE Journal. 1961; 7(2):221-225.

12.   Crane LJ. Flow past a stretching plate. ZAMP.1970; 21645-647.

13.   Raju CSK,.Sandeep N, Sulochana C, Sugunamma V, JayachandraBabu M , Radiation, inclined magnetic field and cross-diffusion effects on flow over a stretching surface. Journal of the Nigerian Mathematical Society. 2015);  34(2):(169-180.

14.   Raju CSK, Sandeep N,. Jayachandra Babu M, Sugunamma V, Dual solutions for three-dimensional MHD flow of a nanofluid over a nonlinearly permeable stretching sheet. Alexandria Eng. J. 2016;55: 151–162. doi:10.1016/j.aej.2015.12.017.

15.   Jayachandra BabuaM, Raju CSK,Sandeep N, SulochanaV, Sugunamma C,. Radiation, inclined magnetic field and cross-diffusion effects on flow over a stretching surface. Journal of the Nigerian Mathematical Society. 2015;34(2):169-80.

16.   SandeepN, Sulochana C, Raju CSK, Jayachandra Babu M,  Sugunamma V.Unsteady boundary layer flow of thermophoretic MHD nanofluid past a stretching sheet with space and time dependent internal heat source/sink.Appl. Math.2015;10(1):312-327

17.   Sohail Nadeem a, Rizwan Ul Haq a, Zafar Hayat Khan,Heat transfer analysis of water-based nanofluid overan exponentially stretching sheet.Alexandria Engineering Journal.2014; 53: 219–224

18.   Reddy JV, Sugunamma V, Sandeep N, Raju CSK.Induced Magneticfield Effect on Stagnation-Point Flow of Magneto-Nanofluids Towards a Stretching Sheet.Advanced Science, Engineering and Medicine.2015; 7(11).968-974

19.   Ishak A, Nazar R, Pop I. Heat transfer over a stretching surface with variable heat flux in micro polar fluids. Phys. Lett. A. 2008; 372(5): 559–61.

20.   Aziz MA, Mixed convection flow of a micropolar fluid from an unsteady stretching surface with viscous dissipation. Journal of the Egyptian Mathematical Society. 2013 ;21:385-394

21.   Raju CSK, Sandeep N, Babu MJ, Reddy JVR. Stagnation-Point Flow of a Ferrofluid Towards a Stretching Sheet.Journal of Nanofluids. 2016;5: 1–8. doi:10.1166/jon.2016.1209.

22.   Mamatha SU, Mahesha, Raju CSK .Cattaneo-Christov on heat and mass transfer of unsteady Eyring Powell dusty nanofluid over sheet with heat and mass flux conditions. Informatics Med. Unlocked. 2017;9:76-85doi:10.1016/j.imu.2017.06.001.

23.   Jumah RY, Fawzi A, Abu-al-rub F, Darcy- Forchheimer mixed convection heat and mass transfer in fluid saturated porous media.  International Journal of Numerical Methods for Heat and Fluid Flow. 2001;11(6): 600 – 618

24.   Hayat T, Haider F, Muhammad T,  Alsaedi A. Three-dimensional rotating flow of carbon nanotubes with Darcy-Forchheimer porous medium.  Plos/one. 2017 ;1–18.

25.   Seddeek MA, Influence of viscous dissipation and thermophoresis on Darcy – Forchheimer mixed convection in a fluid saturated porous media.Journal of Colloid and Interface Science. 2006;293 :137–142. doi:10.1016/j.jcis.2005.06.039.

26.   Animasaun Isaac Lare. Casson fluid flow with variable viscosity and thermal conductivity along exponentially stretching sheet embedded in a thermally stratified medium with exponentially heat generation. Journal of Heat and Mass Transfer Research.2015;2:63-78

 


 

 

 

 

Received on 11.06.2017       Modified on 12.07.2017

Accepted on 18.09.2017      ©A&V Publications All right reserved

Research J. Science and Tech. 2017; 9(3):453-460.

DOI: 10.5958/2349-2988.2017.00079.1